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Abstract. Lazard and Rouillier in [9], by introducing the concept of 
discriminant variety, have described a new and efficient algorithm for 
solving parametric polynomial systems. In this paper we modify this 
algorithm, and we show that with our improvements the output of our 
algorithm is always minimal and it does not need to compute the radical 
of ideals. 
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1 Introduction 

Many problems in science and engineering such as biology, chemistry, computer 
science, robotics, and so on, can be reduced to solving a parametric polynomial 
system. In 1948, Tarski |13j . has published a quantifier elimination method for 
the elementary theory of real closed fields. However the Tarski method provides 
a decision method, which enables one to decide truth any sentence of the theory, 
but the complexity of this approach makes it unsuitable for non trivial problems. 
Later, Boge [7] made significant improvement to this method. We may also cite 
the work of Seidenberg [T^], and later Cohen [5] on this subject. Then Collins 
introduced a completely new method called Cylindrical Algebraic Decomposition 
(CAD), see [3]. 

Weispfenning in [IB] , using the concept of comprehensive Grobner bases, has 
presented an algorithm to describe the complex solutions of a parametric system, 
but his method gives no information on its real roots. Triangular decomposition 
(see [811115114] for example) is another approach that can be used to decompose 
a parametric system. Complexity results for this method can be found in [5]. 

Lazard and Rouillier [S], have introduced the concept of discriminant va- 
riety. Using a projection the parameters space is divided into two parts: the 
discriminant variety and its complement. The discriminant variety is a part for 
which the system has a non-generic behavior (see Section 2 for more details). In 
[TU] . Moroz has proved that the degree of the computed minimal discriminant 
variety is bounded by D := (n — d + £)k^- n ~ d+1 ^ and the variety can be computed 
within the bit complexity c^ c, ( 1 )_D c '( ,l ) if the input system has n — d polynomial 
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equations, I polynomial inequations in n — d variables of degree bounded by k 
with coefficients in a polynomial ring of d parameters with rational coefficients 
of bit-size at most a. 

To compute two component of discriminant variety (critical and singular 
varieties), by the algorithms described in [9], we may need to compute the radical 
and primary decomposition of some ideals. Moreover, the computed discriminant 
variety is not always minimal. In this paper, we give some improvements of this 
algorithm such that to compute these two components, we do not need the 
computation of primary decomposition and radical of ideals. Also, we prove 
that, with our improvements, the computed discriminant variety is minimal (see 
Section 3). 

In Section 2, we recall the definition of the discriminant variety from [3]. 
Section 3 is devoted to computing two components of the discriminant variety 
without using primary decomposition and radical of ideals. In section 4, we 
compare by an example, computing these two components, by the algorithm in 
[9J and our algorithm. 

2 Discriminant Variety 

In this section, we recall the definition of discriminant variety, and we state some 
results from [5] to compute it. For more details see [5J. Throughout this paper, 
we use the following notation. 

Notation 1 Let R — C[tii, . . . , Ud, Xd+i, ■ ■ ■ , x n ] be a polynomial ring where C is 
the field of complex number, U = {ui, . . . , Ud} is the set of parameters and X = 
{xd+i, ■ ■ ■ , x n } is the set of variables. Let also p±, . . . ,p 8 , ft, . . . , ft, £ R be some 
polynomials. Let £ = {pi, . . . ,p s } is the set of equalities and J- — {/i, . . . , fg} is 
the set of inequalities, we define 

C = {x e C" | Pl (x) = ■ ■ ■ = Ps (x) =0,f 1 (x)^0,..., f e (x) ? 0}. 

We denote by Liu '■ C" — > C d the projection map LIjj(U,X) — U . Let LIu(C) be 
the Q-Zariski closure of Liu (C) a nd 5 be its dimension. 

We can easily prove that LIu(C) — LIu(C). 

Definition 1. The affine variety W C LIu(C) C C d is called discriminant vari- 
ety of C w.r.t. Liu if the followings hold: 

1. W^ HuW) . 

2. W = LIu(C) if and only if LI^J 1 n C is infinite for all u £ LLu(C). 

3. LIu{C)\W is union of finite number of connected component of dimension 5 
(i.e. LLuiC) = Uj=i^i where Ui is a connected component of LLu(C)\W and 
dim(W,) = 6 for all 1 < i < k). 

4- LIu{C)\W is a covering space oftii for all 1 < i < k. 

From the fourth property of the above definition, we can conclude that for 
each connected component U C LTu{C)\W, there exist a finite set of indices A 
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and disjoint connected component subsets {VaIag/I of C such that II u l (pl)C\C = 
{JxeA and LIjj : V\ — > U is a diffeomorphism for all AsA 

Now we introduce some other notations that use through this paper. 
Notation 2 

1. Ooo is the set of u £ LIu(C) swc/i i/iai LI^iJA) C~\C is not compact for any 
compact neighborhood hi of u. 

2. O s d is the projection by LIjj of the components of dimension less than 5 of 
C. 

3. O c is the projection by U\j of the critical locus of the union of the components 
of dimension S of C. 

4- Ojr is the projection by H\j of the intersection of C with Yii=i fi = fi- e - 

Or = n u (cnv(uUf i ))). 

5- S i n g is the singular locus of IIu{C). 
6. W x is the Zariski closure of O x . 

The afhne variety Ooo U O se i U O c U Ojr U S i ng is contained in any discriminant 
variety. Now we show that it is the minimal discriminant variety. 

Proposition 1. W D = 0!1 UO S(i UO c U0 7 U O sing = Woo U W sd UW C U%U 
W s i n g is the minimal discriminant variety of C w.r.t. U\j . 

Proof. See [5], Theorem 1. 

Now, the main problem is how we can compute Wd by the above proposition. 
Let 1 c Q[U,X] be an ideal such that V(X) = C. By Pre proces sing (resp. 
PropernessDefects) algorithm in [S] we can compute I, IIu(C), 6 and Wjr 
(resp. Woo)- For computing W st z we may need to compute a primary decompo- 
sition of I, however for the real examples it is empty. In the next section, we 
will propose some improvements of the algorithms described in [S] to compute 
W c and Wsing for computing Wd- 

3 Statement of the main results 

In this section, we state our main results to compute W c and W S ing- For this, 
we need the following notations. 

Notation 3 Let Y = [m, . . . , Ud, %d+i, ■ ■ ■ , x n ] and I C Q[Y] be an ideal. Let 
Y' C Y be a subset of the variables and k < jfY' be a positive integer. We denote 
by Jacy/(I) the ideal generated by minors of size k of the Jacobian matrix w.r.t. 
Y' of a system of generators of L. 

Since Q[U, X] is a Cohen-Macaulay ring, the codimension of I (see the pre- 
vious section) is equal to n — 6. Thus Jacobian criterion and definition of W c 
and Wsing imply that 

W c = V((X + Jac£- 5 (Z)) n Q[U}) 
Wsing - V((ZnQ[(7]) + Jac^ 5 (XnQ[(7])). 
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We will show that to compute W c and W S i„ s , one can replace Jac^- 5 (I) and 
Jac^ZnQ^]) by J ac%~ *({£}) and 3acff\{£) nQ[Z7]) respectively. The com- 
putation of JaCy~ 5 (I) is more expensive than the computation of JaCy _<5 ((£ )) 
for any V C {£/, X}, since the number of component of X is more than the 
number of component of £. 

Lemma 1. The ideals I+JacJ" " <s : (T) and T+JacJ" 6 ((f)) fand (£)+Jac^ _<5 ((£ »j 
ftaue i/ie same zeros outside V(J} i=1 /i). 

Froo/. Let / = ]li=i /<• Since 1 = ( £ ) ■ f°°, then (O C I and therefore 
V(I) C V ((£)). Conversely, let g G I. Then by definition, there exists (eN 
such that gf G (5). Let a G V((£)) and /(a) / 0, so (gf)(a) = 0, then 5 (a) = 
0. This follows that a G V(I) and V((f))\V(/) C V(J). Thus V(2)\V(/) = 
V((£))\V(/). 

Now let I = (<7i, . . . , <?fc). The Jacobian matrix of the ideal I does not depend 
on the chosen generators of I (see [6], Corollary 16.20, page 405). Since (£) C I, 
then Jac£~ ff ((£)) C Jac£ _4 (J), so V(Jac£ -4 (Z)) C V(Jac™-' 5 ((£))). Conversely 
by definition of X, for every gi G {<?i, . . . , <7fe} there exists a natural number 
such that gif li G (£). Set t — max{^ | 1 < i < k}. Then g^f* G (£) for any i, 
and Jacx(/*I) C Jac* ((£)). Thus 

Jw x (fX)=l — ) 1/ . . 

V aX J /fcxd V UX 3 UX 3 J kxd 

Now let a G V(2)\V(/), then &(a) = for 1 < i < fc, so 
Jacx(/*I)(a) = (f(a)^\ 
= /W%^) =/ t (a)Jac x (X)(a) 

thus V(Jac^ 5 ((£)))\V(/) C V(Jac^ 5 (I))\V(/) and it ends the proof. 

The computation of {(£) + J&c x ~ 5 ((£ ))) nQ[J7], in practice, is usually faster 
than the computation of (X + J &c x ~ s ((£))) n Q[U]. So it would be better if 
we replace X by (£) to calculate W c and W s i ng . For instance, in [3], Example 
5.1 (cuspidal manipulators) if we compute ((£) + 3a,c x ~ s ((£))) fl Q[£7] by the 
corresponding algorithm that we have implemented in Maple 12, it takes 11133 
seconds, while for computing (X + Jac x ~ s ((£))) D Q[U], we need 15453 seconds 
(the timings were conducted on a personal computer with Intel(R) Core(Tm)2 
Duo CPU T5670O1.80 GHz and 1.79 GHz, 1.99 GB of RAM). 

The proof of the above theorem implies the following results. 

Corollary 1. The ideals {XnQ[U})+ J ac^ (XnQ[U]) and (XnQ[U])+3acfj- s ((£)n 
Q[U]) have the same zeros outside V((T[Li /*) n ^PY)- 
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Theorem 1. W 00 UW sd UWjUV((l + Ja4^(^)))nQ[[/])UV((lnQ[(/]) + 
J&c^j~ s ((£) DQ[U])) is equal to Wd, o,nd therefore it is the minimal discriminant 
variety. 

It is worth noting that Lazrad and Rouillier in [5], have proved this theorem in 
the case that dim (((I + 3ac x ~ s (T)) n Q [?/])) < 8 (see [5], Proposition 3). Now, 
we describe Critical algorithm to compute W c . 



Critical algorithm 

Input: £, G, 8, U, X where G is the reduced Grdbner basis of I w.r.t. ~<u,x 

where -<u and -<xare degree reverse lexicographical ordering. 
Output: The reduced Grobner basis G c w.r.t. -<u such that Wjr U W c = 

G jac :=The reduced Grobner basis of (G) U Jac^- -<5 ((£ }) w.r.t. -<u,x\ 

G c :=G jac D®[U]; 

Return (G c ); 



We present Singular algorithm to compute W S ing- Note that if 8 = d, then 
W s i n g = to optimize the computation. 



Singular Algorithm 

Input: £ , Gn, 8, U, X where Gn is the reduced Grobner basis of I n Q[U] 
w.r.t. -<u where -<u is degree reverse lexicographical ordering. 

Output: The reduced Grobner basis G S i ng w.r.t. -<u such that WjUW sm9 = 
%UV((G sm5 )). 

if 8 < d then 

G S ing'= The reduced Grobner basis of (Gn) U Jac^ - * 5 ((Gn)) w.r.t. 
else 

G sing {f}: 

end if; 

Return (G slng ); 



4 Example 

In this section by an example we compare the algorithm in [3] and our algorithm 
described in Section III to compute the discriminant variety. In this example we 
will show, whereas by the Critical algorithm in [9] we need to compute radical 
of ideals and the output is not the minimal discriminant variety while by our 
algorithm we obtain the minimal discriminant variety without using radical of 
ideals. 
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Let C — {(a,r,x,y) e C 4 | ax 2 y + bay 3 — r 3 = 0, a - r 2 = 0, r > 0} be a 
constructible set. Thus £ — {ax 2 y + bay 3 — r 3 ,a — r 2 } is the set of equalities, 
J- = {r} is the set of inequalities, U = [r, a] are the parameters and X = [x, y] 
are the variables. By the Preprocessing algorithm we get 

S = 1 

I = (r 2 — a, bay 3 + ax 2 y — ra) 
ll^T) = (r 2 - a) 
Wjr = ar. 

The Output of PropernessDefects algorithm is 

Woo = r 2 - a. 

Using SmallDimension algorithm, we have to compute primary decompo- 
sition of X which is I — (a, r 2 ) n (r 2 — a, x 2 y + by 3 — r) and therefore W s d = 0- 
If we run Critical algorithm in [9], on this constructible set, we get 

G c = r 2 - a 

Gging 

Property = "NeedRadicaF 

Thus, by the algorithm in [9], to calculate the discriminant variety we need to 
compute radical of ideals and the computed discriminant variety is not minimal 
while by our algorithm W c = and W S ing = and the computed discriminant 
variety is minimal. 
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